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ABSTRACT 


A  Coherent  Optical  Subsystem  (COSS)  for  scanning  photographic  Imagery  Is 
Interfaced  to  a  minicomputer  forming  an  Optical  Power  Spectrum  Analyzer  (OPSA) 
designed  for  applications  of  feature  classification  and  "Image  quality" 
measures  on  a  production  basis.  Applications  are  directed  at  finding  non- 
correlation  areas  of  aerial  Imagery  Intended  for  use  In  automated  photo- 
gramnetrlc  compilation  equipment.  The  Image  quality  functions  are  concerned 
with  photogrammetrfc  Information  and  the  overall  utility  of  the  image. 

INTRODUCTION 


The  OPSA  was  designed  after  a  considerable  research  effort  on  the  Recording 
Optical  Spectrum  Analyzer  (ROSA)  had  established  that  optical /digital 
processing  would  be  a  beneficial  means  of  rapidly  screening  aerial  photography 
(Ref.l).  Inherent  In  this  task  was  the  work  of  many  others  who  laid  the 
ground  work,  tested  pattern  classification  algorithms,  and  applied  power 
spectrum  analysis  In  a  variety  of  applications  (Ref.  2,  3,  4,  5). 

Data  acquisition/processing,  system  hardware  design,  and  power  spectrum 
optical  design  have  all  evolved  since  the  Initial  efforts  of  Lendarls  and 
Stanley  (Ref.  5).  Beam  filtering  masks  have  been  replaced  by  a  polar  symmetric 
ring-wedge  detector  which  offers  a  tremendous  advantage  In  data  acquisition 
time.  Two  dimensional  positioning  of  the  Image  with  respect  to  the  Input  beam 
aperture  has  always  been  a  trade-off  between  accuracy  of  the  hardware  and  time. 
Manual  positioning  was  slow,  x-y  positioning  stages  along  with  their  computer 
compatablllty  became  the  next  generation.  Inherent  to  the  stages  were  their 
time  limiting  stabilization  characteristics  and  their  Inability  to  keep  up  with 
the  ever  Increasing  capabilities  of  the  digital  computer. 

Telecentrlc  optical  scanning  systems  represent  the  latest  technology  and 
are  a  true  challenge  to  the  system  architecture  of  any  present  day  computer. 
While  telecentrlc  systems  offer  the  speed.  Its  stringent  design  parameters 
must  Include  optical  stability,  state-of-the-art  optics,  thermal  and  laser 
stability  consideration.  Once  the  optical  design  has  been  optimized,  the 
designer  must  then  confront  the  digital  electronic  design  of  the  scanner  and  Its 
Interface  to  the  computer.  System  applications,  scanning  methodology,  data 
accuracy  requirements,  and  computer-controlled  functions  must  all  be  established 
to  meet  the  needs  of  the  end  user.  Most  Important  Is  optimization  of  the 
electronics  and  Its  direct  relation  to  the  optical  subsystem  and  the  computer. 

Real-time  data  processing  requires  nanosecond  timing,  and  places 
strict  requirements  on  system  architecture*  operating  system  overhead,  special 
Interrupt  handling  with  associated  high-speed  drivers,  and  microprogrammed 
Instruction  subfunctions.  Requirements  Insist  on  a  thorough  knowledge  and 
optimization  of  design  parameters  for  successful  Integration  of  the  above 
topics  Into  a  production  oriented  optical  power  spectrum  analyzer.  The  OPSA 
design  has  Integrated  a  number  of  state-of-the-art  technologies  and  represents 
the  culminated  efforts  of  a  group  of  diversified  professionals. 
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Present  work  on  the  OPSA  has  emphasized  software  design  for  both  testing 
and  production  applications.  Directly  associated  with  this  software  are  the 
high-speed  software  driver  requirements  for  real-time  classification 
processing.  This  paper  will  therefore  address  hardware  functions  of  the 
OPSA  and  software  parameters,  with  a  brief  overview  of  some  preliminary  results. 

APPLICATIONS 

Due  to  the  Importance  of  aerial  photography  to  mapping  and  associated 
stereo  compilation  equipment,  It  was  evident  that  a  screening  device  would 
provide  information  on  the  utility  of  the  Image.  Correlation  of  a  ground 
point  In  a  stereo  pair  of  Images  relies  on  digitally  overlaying  the  points 
via  a  best  fit  of  a  set  of  points,  where  the  set  of  points  Includes  the 
desired  ground  point  (Ref.  7).  Non-correlation  results  In  errors  and  the 
requirement  to  find  another  set  of  Imagery  to  cover  the  non-correlating  area. 
Present  screening  Is  done  manually  to  enclose  possible  adverse  (non-correlating) 
areas  and  In  a  separate  step  to  find  the  percentage  of  cloud  cover  of  an  Image. 

Requirements  for  adverse  areas  are  concerned  with  Identifying  areas  of 
non-correlation  on  stereo  compilation  equipment,  which  Implies  homogeneity  of 
the  area  or  a  density  outside  the  range  of  the  equipment.  Spectral  Information 
(radial  frequency,  angular  orientation)  lends  Itself  directly  to  Identifying 
homogenous  areas,  and  absolute  density  calibration  of  the  OPSA  allows  density 
measurements  to  be  analyzed.  Energy  (Intensity)  distribution  of  the  Fourier 
diffraction  pattern  Itself  Implies  texture  (feature)  characteristics,  with 
homogenous  areas  Implying  little  scattering.  This  Is  a  simplistic  approach 
to  the  problem,  but  serves  as  a  starting  point.  Thus,  there  Is  a  tremendous 
advantage  to  analysis  of  optical  diffraction  patterns  (optical  processing) 
verses  texture  analysis  of  a  set  of  pixel  data  values  generated  by  a  standard 
optical  Imaging  system. 

The  preliminary  adverse  area  feature  set  consists  of  clouds.  Other 
feature  sets  will  consist  of  water,  "dense"  forests,  snow  and  sand.  Image 
quality  parameters  consist  of  density  calibration  of  the  OPSA,  histograms, 
density  range,  contrast  measurements,  and  film  duplication  transfer  functions. 

HARDWARE 

The  COSS  Is  a  computer-controlled  roll  film  scanning  subsystem  which 
images  a  Fourier  diffraction  pattern  on  a  64-element  ring-wedge  detector. 

It  consists  of  three  major  elements:  (a)  a  telecentrlc  scanner  Illumination 
system,  (b)  the  Fourier  transform  optics,  and  (c)  the  precision  continuous 
film  transport.  There  are  six  aperture  sizes  from  Imm  to  10mm,  with  any 
combination  of  1mm  to  10mm  for  across  and  along  film  sample  spacing.  This 
allows  optimization  of  scan  time  per  frame  and  allows  matching  resolution 
requirements  to  meet  specific  applications. 

Very  small  beam  positioning  errors  are  compensated  for  with  a  pair  of 
plezolectrlc  drivers  which  are  driven  via  error  signals  from  a  quadrant  cell, 
which  detects  optical  alignment  errors.  A  galvanometer  provides  a  linear 
scan  across  the  film  plane  and  Is  controlled  via  a  digitally  synthesized  Input 
for  sampling  strategy.  The  wave  form  for  the  galvanometer  Is  stored  on 
prograsiaable  read-only  memory  forming  an  analog  ramp  which  establishes  the 
scan  frequency. 

Synchronization  of  the  data  sample  with  the  scanning  laser  beam  at  the 
film  plane  Is  accomplished  with  an  optical  clock.  A  grating  Is  scanned  via 
the  deflection  of  a  second  laser  beam  off  the  galvanometer.  The  grating 
Interrupts  the  second  laser  beam  at  a  rate  corresponding  to  1mm  Increments  at 
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the  film  plane.  This  Interrupted  beam  Is  focused  on  a  photo  cell  whose  output 
Is  connected  to  logic  circuitry  thus  establishing  synchronization  for  data 
Input  to  the  computer. 

The  desired  aperture  Is  selected  via  a  turret  assembly  of  6  objective 
lenses  In  combination  with  6  beam  masks.  This  assembly  gives  a  diffraction 
limited  spot  of  light  with  a  center- to-edge  Intensity  fall  off  less  than 
50  percent  at  the  film  plane. 

A  64-element  ring-wedge  detector  Is  connected  to  a  High-Speed  Diffraction 
Pattern  Sampler  (HSOPS)  with  each  element  having  Its  own  associated  amplifier 
and  analog  storage  circuitry  In  tire  HSDPS.  An  analog- to- digital  converter 
sequentially  addresses  and  converts  the  analog  signal  every  40  microseconds. 
The  digital  output  represents  the  Intensity  of  the  element  In  a  range  of 
about  32K.  This  digital  output  Is  Interfaced  to  a  minicomputer. 

An  attenuation  wheel  of  8  neutral  density  filters  (0.0  to  1.0)  Is 
Incorporated  to  prevent  saturation  of  the  ring-wedge  detector.  Density 
calibration  of  the  COSS  Is  attained  by  selecting  each  of  the  8  density  filters 
and  correlating  the  Intensity  at  the  detector  with  that  filter's  known 
density.  Access  to  the  attenuated  beams'  Intensity  Is  also  monitored  by  a 
photo  cell  placed  near  the  film  plane.  Another  photo  cell  monitors  the  beam 
before  the  attenuation  wheel,  thus  any  laser  Intensity  drift  during  the 
calibration  procedure,  or  during  the  scanning  mode,  can  be  corrected. 

SOFTWARE 


A  matrix  eigenvector  software  routine  Is  used  as  a  means  of  establishing 
the  merit  of  the  statistics  used  In  feature  classification  decisions.  This 
routine  Is  an  adaptation  of  a  method  published  In  a  paper  by  R.  L.  Mattson 
and  J.  E.  Dammann  (Ref.  8).  It  consists  of  a  technique  to  approximate  the 
eigenvector  corresponding  to  the  largest  eigenvalue  of  a  matrix.  This 
maximized  eigenvalue  matrix  represents  the  weighting  function  (relative 
merit)  of  the  statistics  to  be  used  In  the  classification  algorithm.  Thus  the 
relative  merits  of  a  group  of  statistics  can  be  used  as  a  means  of 
establishing  the  most  useful  statistics. 

An  x-y  digitizer  Is  used  as  a  means  of  establishing  the  "ground  truth" 
of  a  particular  feature  type.  The  digitized  x-y  coordinates  consist  of  a  . 

set  of  contiguous  values  representing  the  circumference  of  a  feature  which  j 

are  passed  through  either  a  three  or  six  dimensional  warping  function.  (Ref.  j 

9,10).  The  Input  for  the  warping  function  consists  of  a  set  of  known  } 

sample  points  Identified  from  the  scanned  data  base,  thus  there  Is  correlation 
(1  to  1  mapping)  of  the  scanned  data  base  with  the  x-y  digitizer  data  base. 

The  correlated  coordinates  are  then  used  as  Input  to  a  software  routine  that 
finds  all  points  within  the  circumference  enclosed.  The  software  Is  flexible 
enough  to  allow  Input  of  a  set  of  circumferences  or  lines  (non-enclosed).  This 
"around  truth"  Is  then  compared  (overlayed)  to  the  output  of  a  feature 
classification  algorithm  used  on  an  OPSA  scanned  data  base.  Thus,  the  percentac 
overlay  of  the  two  data  bases  gives  a  relative  accuracy  of  the  decision 
algorithm. 

The  x-y  digitizer  Is  also  used  as  a  means  of  Inputlng  D-log  E  curves  for 
generating  transfer  functions  used  for  duplication.  The  transfer  function  Is 
generated  from  the  D-log  E  curve  of  the  original  Image  and  the  D-log  E  curve 
of  the  proposed  duplicating  film.  It  Is  assumed  that  some  of  the  Image 
"characteristics"  (density  range,  contrast,  histogram)  of  an  original  Image 
can  be  representative  of  the  Image.  By  passing  the  OPSA  Image  data  through 
the  transfer  function,  the  "characteristics”  of  the  duplicated  Image  can  be 
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made,  thus  determining  If  the  duplicate  meets  Its  specific  requirements.  In 
a  similar  manner.  If  the  desired  "characteristics”  are  known  then  all  D-log  E 
curves  In  the  data  library  can  be  searched  to  find  a  best  fit,  giving  the 
duplicate  film  type  and  processing  requirements. 

Computer  control  of  the  COSS  allows  the  setting  of  the  along  and  across 
film  sample  spacing,  galvanometer  position,  computer  or  manual  mode,  film 
width,  and  slew  direction  and  speed.  The  frame  counter,  film  length  counter, 
laser  power,  optical  clock  status  can  all  be  monitored  via  the  three  software 
drivers  that  control  the  OPSA. 

RESULTS 

Classification  results  are  based  on  relating  the  x-y  digitizer  ground 
truth  to  the  decision  classification  algorithm  results  for  cloud  features. 
This  allows  a  total  percentage  comparison  and  a  point  by  point  assessment 
of  the  accuracy  of  screening  for  clouds.  A  two-dimensional  feature  vector 
consisting  of  a  threshold  of  the  D.C.  Intensity  and  a  ratio  of  a  sum  of 
low-frequency  spectral  rings  to  higher- frequency  spectral  rings  were  used 
for  classification.  Selection  of  the  threshold  and  sum  of  ring  elements 
was  determined  by  examining  ring  combination  ratios  and  using  the  combination 
with  the  greatest  relative  separation.  Tests  were  run  under  varying  condi¬ 
tions  for  2,  4,  and  6  mm  apertures  of  a  test  Image.  Preliminary  results 
Indicate  better  than  90  percent  accuracy. 
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